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When an electrically charged system is subjected to the action of an electromagnetic field, it 
responds by generating an electrical current. In the case of a multicomponent plasma other effects, 
the so called cross effects, influence the flow of charge as well as the heat flow. In this paper 
we discuss these effects and their corresponding transport coefficients in a fully ionized plasma 
using Boltzmann's equation. Applications to non-confined plasmas, specially to those prevailing 
in astrophysical systems are highlighted. Also, a detailed comparison is given with other available 
results. 



I. INTRODUCTION 

In this paper we address ourselves to the study of the behavior of a fully ionized dilute plasma 
which is subjected to the action of an electromagnetic field. Although the primary response of 
the plasma to an electric field is the generation of an electrical current, if a magnetic field is also 
present other effects arise including a heat current. This gives rise, according to the tenets of 
classical irreversible thermodynamics (CIT) to a number of interesting and important cross effects 
which, to our knowledge, have received little or no attention in the literature when examined under 
the framework of microscopic equations. This is the main objective of our work, to use the well 
known Boltzmann equation to calculate the explicit forms of the flux-force relations demanded 
by CIT and provide explicit expressions for the ensuing transport coefficients as functions of the 
density, temperature and magnetic field. We believe that the results will turn out to be valuable 
for many astrophysical systems and also other situations related to non-confined plasmas. 

To achieve this program, in Section II we summarize the kinetic basics of the problem using 
standard techniques of the kinetic theory of gases, in Section III we establish all flux-force relations 
appearing in the system including the explicit calculation of the transport coefficients, and in Section 
IV we present a critical discussion of our results including a comparison with others available in 
the literature. 
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II. KINETIC BACKGROUND 

Foundations of the methods here applied are available in the somewhat extensive literature on 
the subject 0, 0, H 0, 0, H so we shall avoid detailed algebraic steps which the reader may 
easily find in any of the works here cited. 

We consider a dilute binary mixture of charged particles with masses m a , rrib; charges e a — 
— eb = e taking Z = 1 merely for didactic purposes. The number density n = n a + rib and the mass 
density p = m a n a + rribrib where n$ is the number density of species i (i = a, b). The fields E and 
B acting on the system are the self-consistent fields as determined from Maxwell's equations and 
the possibility of having an additional external magnetic field is taken care off by defining the total 
field B T = 5 + 5 (ext) . These fields are included in Boltzmann's equation through the Lorentz force 
j§] so that if fi (r , Vi , t) is the single particle distribution function for species % this equation reads 
as 

Here, the electric field is included as part of the total conservative force. 

Pi = ^ (ext) + eS (2) 

and 

J ififj) = J- J {ft (^) fj fij) - fi fj ° (vi, -> ^, tfj) dVj dv[ dVj (3) 

where the primes denote the velocities after a collision (Vi, Vj — » v' i , Vj) and the probability that 
this collision occurs satisfies microscopic reversibility namely, 

a (vi, Vj —> v 1 ^ Vj) = a Vj -> v i} Vj) i, j = a, b (4) 

In words, Eq. Q guarantees the existence of inverse collisions among the two species. 

For our purposes, we concentrate on a general feature of Eq. ([1]), namely the derivation of the 
conservation equations for the chosen local variables: pi (r, t) the local mass density for species 
i, pu(r, t) the local baricentric momentum per unit volume and the local internal energy density 
e (r, t). In fact multiplication of Eq. {1J by m,, rriiVi and ^iriivf and integration over Vi leads 
effortlessly to the following results: 

i) Mass conservation equation, 

^ + V • ( Pi u) = -V-Ji (5) 
where u, the barycentric velocity, is defined as 

b 

pu(f, t) = PA (r, t) (6) 

i=a 

with 

u% (r, t) = — j fi (r, Vi, t) Vidvt = (vi) (7) 
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and clearly 



{4>) = — I fi (7, Vi, t)ip (vi) dvi (8) 



11, 



The diffusive flux Jj is 

and the thermal or random velocity Cj is simply given by c, = Mj — u (f , i) . From these definitions 
it follows readily that J a + Jb — Pa (c a ) + pf, (q,) = 0. 
ii) Momentum conservation equation, 

(pit) + V ■ (r fe + puu) = rnFi + Jt x B (10) 

i 

where the total charge current Jt is defined as 

b 

Jt = Yn l e i (v l ) (11) 

i=a 

and the stress tensor is given by 

b 

r k = ^ m i / /' *) ^Qdwi (12) 



Notice that according to Eq. 1(12)1 r fc is a symmetric tensor. Also, in Eq. (fit))) Fj is the total conser- 
vative force acting on species i including the electric field E. 
iil): Internal energy equation, 

p^- t (^j +\7 ■ J Q + r k :X7u- J c - E' = (13) 

where 

pe{f, t)=Y,\pi(4) (14) 

i 

The heat flux Jq is given by 

and the conductive current J c is defined as 

J c = 2_,n i ei{c i ) (16) 
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Finally, E 1 = E + u x B is the "effective" electric field that is, the field measured by an observer 
moving with the baricentric velocity u. 

It is important to underline the fact that the set of conservation equations © , flO|) and (fl3|) is 
incomplete. Indeed we have six independent variables and six equations, but all the currents are 
unknown, Jj, J c , Jq and r fe . In order to derive a complete set of equations we need to express these 
currents as functions of the local state variables which implies clearly that we must solve Eq. (TJ 
for fi (r, Vi, t) and somehow introduce the local variables in the solution. 

But there is another ingredient that we must require form our solutions. We want to write the 
two important currents namely Jq and the conduction current J c in terms of the gradients in the 
system. This implies seeking relations such as 

J Q = -kVT-TV^ 

J c = -6VT-crV0 (17) 

since E = — V0. These are the linear flux-force relations as demanded by CIT where the transport 
coefficients k, a, T and B have to be determined. Here k and a are the ordinary thermal and 
electrical conductivities whereas T and B are the so called Thomson's and Benedicks coefficients 
respectively. Their form will be radically modified by the presence of a magnetic field contrary to 
the case where B = and only an electrical field is present. Notice also that from the definition of 
J , the fact that J a + Jb = and the definition of J c , 

j m a + m b 7 

Jc = e a J a (18) 

m a mb 

so the mass and conduction currents are proportional among each other. We shall not pursue 
diffusive effects here as they have been reported elsewhere [10] . 



III. THE TRANSPORT COEFFICIENTS AND LINEAR LAWS 

To achieve the program outlined in the previous section, we proceed with the solution of Boltz- 
mann's equation. The first stage of this task is standard in the literature [3,0,0,0,0,0,0,0] so we 
merely outline the main steps. One assumes that the single particle distribution function may be 
expanded in a power series of the Knudsen parameter e which is a measure of the spatial gradients 
of the local variables in the system. This series is taken using the local Maxwellian distribution 
as a reference state since this function is, as well known, a solution to the homogeneous term, the 
collisional term in Eq. (TJ. Moreover this expansion holds for times much longer than the collision 
times so that the distribution function is assumed to be a time independent functional of the con- 
served densities. Known in the literature as the Hilbert-Chapman-Enskog expansion it thus reads 
as 

fi (r, vl, t) = f[ 0) (r, [l + e^ (1) (r, vl\) + ...] , (19) 

here is the local Maxwellian distribution function and the dash in all functions in the right side 
of Eq. lfT9|) implies the time dependence through rii (r, t), u(f, t) and e (r, t). 

After Eq. lfl9|) is substituted back into Eq. ^ and recalling that the zeroth order in e term is just 

Sj=a J (fifj) — 0, whose solution is f- °\ one obtains the linearized Boltzmann equations for ip"p 
whose solutions are known to be of the form 



f ] = Mi : Vu + &i ■ V In T + • d io i, j = a, b (20) 
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Equations i|20p . the first order in e solutions to Eq. (TJ characterize the well known Navier-Stokes- 
Fourier regime of magnetohydrodynamics. In Eq. (|20| 



dij = V 



fUi\ n.;nj ( m-i — mi\ Vp rt,n,- , , ~ , 

— + - 1 — + (m iej - rrijei) E' = -djt , 21) 

\n J n \ p ) p pp 



where p = nkT is the local pressure and the temperature is introduced through the standard ideal 
gas relationship, namely, e (r, t) = |fcT(r, t). The vector is called the diffusive force and 
contains two main contributions, the first two terms related to ordinary and pressure diffusion and 
the third and fourth terms, the diffusive effect arising from the of electromagnetic fields. Also, in this 
expression we have assumed that no external forces are acting on the plasmas, i^ ext ) = Q, j = a, b. 

In what follows we shall ignore the tensorial term Bj : Vu since by Curie's principle [ll(, it will 
not couple with first rank tensors such as VT and E', Also, the structure of the vectors A, and 
Dj and is not arbitrary. They must be linear combinations of the independent vectors available, Ci, 
Ci x B and fcj x B^j x B. Therefore, 

Ai = A^&i + Af ) c l x B + Af ) (&iX Bj x B , (22) 

and a similar expression for Dj. The scalar functions A^ 1 "*, A^ and A^may depend on all scalar 

quantities such as cf, B 2 7 (ci ■ BJ , i%i,T and so on. The next step is to expand these functions in 
terms of a complete set of functions namely, the Sonine polynomials. For instance, 

oo 

A (D = £ a WM 5 M (c 2) f (23) 

and similar expressions for A 2 l , etc. Use of Eqs. ([22| and (f23]) in Eq. |20l) allows a formal computa- 
tion of the fluxes in the system. Let us assume that nib >• m a , nib and m a associated with the mass 
of a proton and an electron, respectively, in a hydrogen plasma. Since then, J c — ^-J a we find, 
using Eqs. ((9|), (f20|) . (|22|) . and (f23]) as well as the orthogonality conditions of Sonine polynomials 
that 

m a I 1 v 

+T 



+ 4 2)(0) ^ x B + (S • Bd^] } , (24) 



This result is quite interesting. Out of the infinite number of coefficients required to characterize the 
scalar functions Aa , (j = 1, 2, 3) as illustrated by Eq. (|23|) . only one coefficient is necessary 
for each of these functions in order to compute the conduction current. These coefficients can 
be obtained from the solution of the integral equations that such functions must obey and whose 
solutions will not be discussed here. The procedure to obtain them is, once more, a standard 
technique in kinetic theory and we refer the reader to the available literature if he wishes to purpose 
the details [3, 0,0, 0,(1,0,0, 13] ■ In Eqs. {24]) is the electromagnetic contribution of the diffusive 
force 

i{e) _ n a n b (m + m b ) e ~, 

d * b ~ Jp E (25) 
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To obtain the expression for the heat flux we must first recall a result from CIT. For a multicom- 
ponent system the standard definition of heat flux is (lll | 

JQ = JQ-J2 h ^ ( 26 ) 

r r 

where h r is the enthalpy of the r th component which is equal to \kT for an ideal gas. In kinetic 
theory the analog of this equation is clearly given by 

i—a ^ ' 

where use has been made of Eqs. ((9|), lfl5|) . and lfl9|) recalling that J'q = in the local equilibrium 
state. 

Now, when Eq. (|20| is substituted in Eq. (|27| the term Aj • VT gives rise to the ordinary Fourier 
heat conduction. With all its modalities arising from the presence of a magnetic field, this has been 
thoroughly studied in a separate paper [l2j so we shall not deal with it here. However the term 
from Eq. lf25|) of our interest here, is the one giving rise to the electromagnetic influence on the heat 
current and reads as, 

4 E) = \ m 2 f £ + tf m B x «g + (B • dg) B) (28) 



To analyze Eqs. 11241) and f28j) we consider a cartesian coordinate system and let B point along the 
z-axis, B = Bk, k being the unit vector. Then for any scalar N, we define, 

V||7V=— ; V ± N=— 1 + — j ; V s iV = — j - — i 29 
11 oz ox oy ox oy 

the last definition corresponding to a vector perpendicular to both the parallel and perpendicular 
components, respectively. Let us now apply Eq. lj29j) to our currents. The thermal contribution to 
the conduction current reads 

J(T) = ^|^(l)(0) +fl 2 a (3)(0)j V||r + a (l)(0) VxT + Ba (2)(0) VsT | (30) 

This is the first important result in this paper, since it contains all the thermoelectrical effects in 
the plasma when a magnetic field is present. Notice that if B = 

4 T) =-TVT (31) 

as required by CIT. T is the standard Thompson thermoelectric coefficient which turns out to be 
given by 

T = 5.334^r (32) 

where, t, the mean collision time calculated from the collision integrals is shown to be given 

by 

_ 247r 3 / 2 eg v ^ 3/2 



mpe 4 



(fcT) J/z (33) 
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and eo = 8.554 x 10 12 F/m is the electrical constant. In the presence of a magnetic field there 
are, besides the current parallel to the z-axis unaffected by the field, two additional contributions. 
Firstly a current along Vj_T, perpendicular to B, with a "perpendicular Thomson coefficient" 



n n keT 



where x = ujt, oj — f^-, Larmor's frequency, r is given by Eq. f|33[) , tp is the so-called logarithmic 



m a Ai (x] 

= UJT i UJ = Larm( 

function 

ip = In 

and Ad is Debye's length defined by 



r 3.i x lo-v 



68x 2 + 8.53) 



(34) 



\ ~ 2 



Ad 



A, 



\ ne J 



1/2 



In Eq. (|34|) we also used the result that 

„(l)(o) _ T 



3.1 x lCTV + 68x* + 8.53) 



Ai(x) 



with 



Ai (a;) = 4.7 x 10" 4 a; 6 + 100x 4 + 32x 2 + 1.6 



(35) 



(36) 



(37) 



(38) 



as obtained from the solution of the integral equations [8|, [lOj . 

Finally, there is a third current in a direction both perpendicular to B and to Vj_T which 
resembles a "thermal Hall effect". The corresponding transport coefficient is 



with 



Ts = !^!i Ba mo) 

m a 

Ai (x) v ' 



(39) 



(40) 



To get a clear idea of how these effects compare to each other, the three coefficients are plotted in 
Fig. 1 as functions of x. Emphasis must be made on the fact that a; is a function of n, B and T so 
that it must be handled with care when seeking orders of magnitude. 

Next we examine the influence of the effective electromagnetic field on the conduction current. 
Using the decomposition given in Eq. f29|) one may write that 



J, 



(e) _ 



e n a kT 



{ (4 1)(0) + s 2 rfi 3)(0) ) T Ez k + 4 1)(0) (d£) ± + M 2)(0) (4 e) )J ( 41 ) 



where 



■tb 



(E x +u y B)i + (E y -u x B)j 
(E x + u y B) j - (E y - u x B) i 



(42) 
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and 



r 



n a nb (m a + m&) e 



VP 



(43) 



Using now the fact that mj, <C m a = m e , p = nkT, the definition of p and assuming a fully 
ionized plasma so that n a = rib — §, Eq. (f43~)) may be rewritten as, 



=a n E z k + a ± {d$ 



(44) 



which is the second important result in this paper. First notice that if B = 0, from Eq. (|44j) . 
(T|| = <tj_ = c and since i? = — V<^> 



= - CT ||V0 

This is the well known form for Ohm's law where 

eTn a kT rf(1)(0) 



II 



(45) 



(46) 



which, under the conditions mentioned above and using the value of di 1 ^ "* = 4.31r extracted from 
the solution to the integral equations, gives 



25.92,3/^ 



/m^e 2 ip 



(47) 



This is the ordinary electrical conductivity as, aside from minor numerical differences, first derived 
by Spitzer fljjl US . [l5| and later by Braginski [16] and Balescu [§]. The transversal electrical 

conductivity, when the value for is used, turns out to be 



6n 3 / 2 e 2 (kT) 3/2 /9.52 x 10 4 a; 4 + 20.6a; 2 - 13 



/m a e 2 rp 



A 2 (x) 



and cr s , which corresponds to a "Hall like effect", turns out to be 

_ 67r 3 / 2 eg (kT) 3/2 ^0.109a; 5 + 2.18 x 10 4 a; 3 + 965a; 



/m a e 2 tp 



A 2 (x) 



where 



A 2 (x) = 1.61a; 6 + 3.5 x 10 5 x 4 + 1.65 x 10V 



3.01 



(48) 



(49) 



(50) 



To appreciate the order of magnitude of these coefficients they are plotted in Fig. 2 as functions 
of x. Moreover it is worth pointing out that from Eq. (|48|) it is easily verified that lim^^o °± = cr]| 
and linis^o &s = 0, as seen from Eq. (J49j) . 

The last item we want to deal here with concerns the influence of E' on the heat current which, 
according to Eqs. ([28]), ([29]) and ([42]), is given by 



j7 } = I (*D a E ^ + ^X (1)(3) ) e z % + («r«) ± + *4 2)(1) (4 e) ) J (3D 
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If B — and assuming a fully ionized hydrogen plasma, using the value for obtained from 

the solution to the integral equations one finds that 

JJ) = -B,V0 (52) 



(^f (53) 



where 



e 2 ip \ m e/ 

the so-called electropyrosis or Benedicks effect [T3|. When B ^ and using the results for the 
coefficients ^ and o?- 2 ^ 1 ^ one obtains for the two other coefficients, namely 

157r 3 / 2 e£m e fkT\ 5/2 5.04 x 10~ 3 a; 4 + 1.09 x 10 3 a; 2 +3.2 , s 

B ^= W WJ M^) (54) 



157r 3 / 2 e 2 m e fkT\ 5/2 0.0295a; 3 + 222a; 



and YmiB^o B± = B\\ and also 

e^-i/> \m e J A 2 (ar) 

A2 (a;) defined in Eq. (|50|l . This effect has been systematically ignored in the study of transport 
processes in plasmas. The importance of the three coefficients is exhibited in Fig. 3. 



IV. DISCUSSION OF THE RESULTS 



The results obtained for the nine transport coefficients derived in the previous section is important 
for several reasons, mainly because with one single exception, they have been ignored in the study 
of transport processes in plasmas. Let us begin with the thermoelectric or Thomson coefficients. 
To our knowledge the only other derivation of these quantities is that given by Balescu in his 
excellent monograph in the subject Nevertheless his derivation starts from Landau's version 
of Boltzmann's equation which is a Fokker-Planck type equation This implies in words that 
transport processes in plasmas may be visualized as diffusive processes. Moreover the solution of 
that kinetic equation is carried out using a moment-like expansion a la Grad which is different in 
spirit from the Hilbert-Chapman-Enskog expansion []j|[2(3]. Here we use the latter to solve the full 
Boltzmann equation, which in turn is not a diffusive like approximation. In Fig. 4 we compare our 
results for T± and T s with those of Balescu showing explicitly that although qualitatively similar 
there are differences in the two methods. 

For the electrical conductivit y th e story is quite different. As mentioned earlier, it was first 
derived by Spitzer et al in 1950 [l3j] using a diffusion equation and later by Spitzer and Harm [14j 
from a Fokker-Planck equation similar to the one used in his previous paper, but making a mild 
attempt to obtain the linear-constitutive laws. In both cases they obtain, for B = 0, the (fcT) 3 ^ 2 
dependence as in Eq. (|47| with minor numerical differences in the respective coefficients. Later, 
Balescu repeated this calculation methodically [9j| and obtained the three conductivities for the 
case. 

Finally, the Benedicks coefficients are completely new, to our knowledge they have never been 
reported elsewhere. This leads us to the final remark in this work, namely, when formulating 
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magnetohydrodynamic equations it is convenient to assess if, in the range of densities, tempera- 
tures and values of the magnetic fields, these cross effects may be neglected. Otherwise important 
consequences will be absent in the ensuing results. 
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Figure Captions 

Fig. 1. The values of the three thermoelectric coefficients for n = 10 21 cto -3 and T = 10 7 K . 
The solid line is 7y , the dashed line T± and the dotted line T s . 

Fig. 2. The values of the three electrical conductivities for n = 10 21 cto -3 and T = 10 7 K. The 
solid line is an, the dashed line a± and the dotted line a s . Notice the predominance of 
er s over cry for a certain range of values of x. 

Fig. 3. The Benedicks coefficients for n = 10 21 cm~ 3 and T = \0 7 K. The solid line is B\\, the 
dashed line B± and the dotted line B s . 



Fig. 4. 



The ratio T±/T\\ (dashed) and Ts/T\\ (dotted) is plotted (right) and compared with the 
corresponding Balescu's results (left), as in Ref. 
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